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ABSTRACT 

The low-redshift evolution of the intergalactic medium is investigated using hydrody- 
namic cosmological simulations. The assumed cosmological model is a critical density 
cold dark matter universe. The imposed uniform background of ionizing radiation has 
the amplitude, shape and redshift evolution as computed from the observed quasar 
luminosity function by Haardt & Madau. We have analysed simulated Lyman-a: spec- 
tra using Voigt-profile fitting, mimicking the procedure with which quasar spectra are 
analysed. Our simulations reproduce the observed evolution of the number of Lyman- 
a absorption lines over the whole observed interval of z = 0.5 to z = 4. In particular, 
our simulations show that the decrease in the rate of evolution of Lyman-a absorption 
lines at z < 2, as observed by Hubble Space Telescope, can be explained by the steep 
decline in the photo-ionizing background resulting from the rapid decline in quasar 
numbers at low redshift. 
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1 INTRODUCTION 

Neutral hydrogen in the intergalactic medium produces a 
forest of Lyman-a absorption lines blueward of the Lyman- 
a emission line in quasar spectra (Lynds 1971). Observations 
of quasars at redshifts z > 2 show that there is strong cosmo- 
logical evolution in the number of Lyman-a lines, which can 
be characterised by a power law dN/dz oc (1 + z) 1 (Sargent 
et al. 1980), where iV is the number of lines above a thresh- 
old rest-frame equivalent width W (typically W > 0.32A). 
Studies at high resolution using the Keck telescope find 
7 = 2.78 ± 0.71 for 2 < z < 3.5 (e.g. Kim et al. 1997). 
At even higher z the evolution is still stronger: Williger et 
al. (1994) find 7 > 4 for z > 4 using the CTIO 4m telescope. 
In contrast at low redshifts, observations using the Hubble 
Space Telescope find much less evolution, 7 = 0.48 ±0.62 for 
z < 1 (Morris et al. 1991, Bahcall et al. 1991, 1993, Impey 
et al. 1996). 

Recently, hydrodynamic simulations of hierarchical 
structure formation in a cold dark matter (CDM) domi- 
nated universe have been shown to be remarkably success- 
ful in reproducing this Lyman-a forest in the redshift range 
z = 4 — » 2 (Cen et al. 1994, Zhang, Anninos & Norman 1995, 
Miralda-Escude et al. 1996, Hernquist et al. 1996, Wadsley 
& Bond 1996, Zhang et al. 1997, Theuns et al. 1998). These 
simulations show that the weaker Lyman-a lines (neutral 
hydrogen column density N^{\ < 10 14 cm -2 ) are predomi- 
nantly produced in the filamentary and sheet-like structures 
that form naturally in this structure formation scenario. Ve- 



locity structure in these lines is often due to residual Hubble 
flow since many of the absorbing structures are expanding. 
In contrast, the stronger lines (ATjjj > 10 16 cm -2 ) tend to 
occur when the line of sight passes near a dense virialised 
halo. 

Over the redshift range investigated in these simula- 
tions a photo-ionizing background close to that inferred from 
quasars (Haardt & Madau 1996) is required to explain the 
properties of the Lyman-a forest. In fact, although most 
simulations have assumed a critical density, scale-invariant 
CDM universe, other variants of the CDM model provide 
acceptable fits with relatively small changes to the ionizing 
background (Cen et al. 1994, Miralda-Escude et al. 1996). 
The general success of CDM-like models in explaining the 
high redshift (z > 2) properties of the Lyman-a forest is 
impressive. 

In this Letter we investigate using hydrodynamic simu- 
lations whether a CDM universe with a photo-ionizing back- 
ground dominated by quasars can explain the observed tran- 
sition in the cosmological evolution of the number of Lyman- 
a lines at 2 < 2. 



2 SIMULATION 

We model the evolution of a periodic, cubical region of a 
critical density Einstein-de Sitter universe (Q — 1, = 0). 
We use a simulation code based on a hierarchical P3M im- 
plementation (Couchman 1991) for gravity and smoothed 
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particle hydrodynamics (Lucy 1977, Gingold & Monaghan 
1977, see e.g. Monaghan 1992 for a review) for hydrodynam- 
ics (Theuns et al. 1998). The comoving size of the simula- 
tion box is L/(2h) Mpc, where the Hubble constant today 
is written as Ho = 100ft km s _1 Mpc -1 . We will assume 
h = 0.5 throughout and describe simulations with L = 5.5 
and 22.22Mpc. A fraction tt B h 2 = 0.0125 of the matter den- 
sity is assumed to be baryonic, consistent with limits from 
nucleo-synthesis (Walker et al. 1991, but note the continu- 
ing debate on the deuterium abundance derived from quasar 
spectra which might favour higher values Q B h 2 = 0.019, see 
e.g. Buries & Tytler 1997 and references therein). The rest 
of the matter is assumed to be in the form of cold dark 
matter. To generate initial conditions for the simulations at 
the starting redshift z — 50 we use the fit to the cold dark 
matter linear transfer function from Bardeen et al. (1986) 
and normalise it such that the linearly extrapolated value 
of as — 0.7 at the present day. These simulations use 64 3 
particles of either species hence the gas mass resolutions are 
1.5xl0 8 M Q and2.2xlO 6 M for L = 22.22 and L = 5.5Mpc, 
respectively. 

Gas in these simulations is ionized and photo-heated 
by an imposed uniform background of ionizing photons as- 
sumed to originate from quasars as computed by Haardt 
& Madau (1996). This flux is redshift dependent, mimicking 
the evolution of the quasar luminosity function. Gas can also 
cool by interacting with microwave background photons and 
through collisional cooling. The detailed rates for all these 
processes as a function of temperature are taken from Cen 
(1992) with some minor modifications. We assume ioniza- 
tion equilibrium throughout and use a helium abundance of 
Y — 0.24 by mass. We have compared in detail the results 
from our code against the published results of TREESPH 
(Hernquist et al. 1996, Croft et al. 1997) and find excellent 
agreement (see Theuns et al. 1998 for more details of these 
comparisons and for a description of our code and cooling 
rates). 

The effective mean optical depth f from the simula- 
tions, with this set of parameters, is significantly lower than 
the observed value. Consequently, we have reduced the am- 
plitude rxjj of the ionizing radiation given by Haardt & 
Madau (1996) by a factor of two. Here, Fjjq is the ampli- 
tude of the radiation spectrum at the hydrogen Lyman edge. 
Since f scales approximately as r cx {n B h 2 ) 2 /hF m (Rauch 
et al. 1997), we would obtain the same results by keeping the 
original value of Tyj from Haardt & Madau but increase 
Q B h 2 from 0.0125 to 0.0177, which is still well within the 
range allowed by nucleo-synthesis. 

At several output times we compute simulated spectra 
along lines of sight through the simulation box. Each spec- 
trum is convolved with a Gaussian with FWHM = 8 km s , 
then re-sampled onto pixels of width 3 km s _1 to mimic the 
instrumental profile and characteristics of the HIRES spec- 
trograph on the Keck telescope. Artificial noise is introduced 
by adding a Gaussian random signal with zero mean, and 
standard deviation a = 0.02 to every pixel (i.e. a SNR of 50 
for pixels at the continuum). The absorption lines in these 
mock observations are then fitted with Voigt profiles using 
an automated version of VPFIT (Carswell et al. 1987). 
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Figure 1. Simulated spectra at redshifts 3,2,1 & 0.5 (top to bot- 
tom). We have combined eight spectra through the simulation 
box to extend the wavelength coverage. The chosen sight-lines 
are the same at each redshift. The wavelength range is chosen so 
as to show the same comoving region at every redshift. 



3 RESULTS AND DISCUSSION 

We show in figure |l] examples of simulated spectra at z = 3, 
2, 1 & 0.5 for the L = 22.22 Mpc lower resolution simulation. 
Fluctuations in neutral hydrogen density, caused by gas trac- 
ing dark matter potential wells, produce absorption features 
similar to those in observed spectra. At low redshifts, there 
are large regions of the spectrum with very low absorption. 
These regions are separated by prominent absorption fea- 
tures, most of which are just a single strong line. At higher 
redshifts, there is considerable absorption over most of the 
spectrum and many lines are blended. Many of the strong 
lines at z — 0.5 can be traced back to higher redshifts. There 
is a clear decrease in the comoving number of lines with de- 
creasing redshift. 

The evolution of the column density distribution with 
redshift is illustrated in figure |^. There is clear evolution 
in the simulated column density distribution with redshift. 
The rate of change depends on column density, with higher 
column density lines undergoing stronger evolution leading 
to steepening of the distribution. At z — 2, the column den- 
sity distribution is cx -Wjjj whereas this has steepened to 
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Figure 3. Evolution of the number of lines with given range in column density with redshift from simulations compared against observed 
evolution. Column density cut 10 131 cm~ 2 < A^pjj < 10 14 cm -2 - simulations: circles connected with dotted line; data: open triangles 
(Kim et al. 1997). Column density cut 10 13 ' 77 cm" 2 < N m < 10 16 cm - simulations: circles connected with solid line; data: filled 
squares (Bahcall et al. 1993), open squares (Impey et al. 1996), filled triangles (Kim et al. 1997), filled pentagon (Lu et al. 1996). Column 
density cut 10 14,5 cm" 2 < -/Vjjj < 10 ie cm" 2 - simulation: circles connected with long dashed line; data: long dashed line shows 
evolution from Williger et al. (1994). Large filled circles are simulation results for the low resolution, large box size, run, open circles 
are for a higher resolution, smaller box size, run. Large open pentagon: re-analysis of simulation at z = 0.5, but imposing the ionizing 
background appropriate to z = 2. Data were taken from figure 2 in Kim et al. (1997), except for the Wiliger et al. (1994) data. 



oc A^ 2 - 1 at z — 0.5 (see figure ^). The rate of evolution 
also depends on redshift, with considerably stronger evolu- 
tion at higher redshifts. The number of weak Lyman-a lines 
(< 10 13 ' 1 cm -2 ) remains approximately constant. At red- 
shifts 3 and 2, there is good agreement between the simu- 
lated column density distribution for our higher resolution 
simulation (L = 5.5 Mpc) and the observed one. 

The drop in the number of lines with redshift can be 
quantified by counting the number of lines within a given 
column density range. The simulation results are compared 
to observations in figure []. The simulations reproduce well 
the number of lines at a given redshift for all three column 
density cuts. Crucially, they also match very well the ob- 
served number of lines at low redshift. Consequently, the 
hierarchical picture of galaxy formation in a critical density 
universe, coupled with the observed evolution in the quasar 
luminosity function, can explain the observed evolution of 
the number of Lyman-a lines over the entire observed red- 
shift range. 

We have re-analysed the z = 0.5 output time after in- 
creasing the imposed ionization flux from its z = 0.5 value 



to the value appropriate to z — 2. The number of lines with 
10 13 ' 77 cm" 2 < JVttt < 10 16 cm" 2 is shown by the open 
pentagon in figure H. This point falls onto the extrapolation 
for the number density evolution for z > 2. Consequently, 
the dominant reason for the higher number of lines at low z 
compared to what would be expected by extrapolation from 
high z, is the decrease in ionizing flux from z = 2 to z = 0, 
itself a consequence of the evolution of the quasar luminosity 
function. 

An estimate of the reliability of these simulations can be 
obtained by comparing the two simulations run at different 
resolutions. The higher resolution simulation produces more 
lines at all column densities, but the difference between the 
two simulations is well within the error bars of the observa- 
tional results. This gives us confidence that we can reliable 
predict the number of lines from these simulations. 

In summary: our numerical simulations show that the 
properties of the Lyman-a forest are in excellent agreement 
with what is expected in a cold dark matter universe with 
a photo-ionizing background dominated by quasar light. In 
particular, our simulations show that the observed decrease 
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Figure 2. Evolution of the column density distribution with red- 
shift, for z = 3, 2, 1, 1.5 and 0.5. Distributions for 2 = 2 and 
2 = 3 are from the higher resolution (L = 5.5 Mpc) simulation, 
the others are for the lower resolution (L = 22.22 Mpc) simula- 
tion. Observational data from Hu et al. (1995, HKCSR), Petitjean 
et al. (1993, PWRCL) and Kim et al. (1997, KHCS) are shown 
for comparison, and their mean redshifts are indicated. Scalings 
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in the rate of evolution of Lyman-a absorption lines at z < 2 
can be explained by the steep decline in the photo-ionizing 
background resulting from the rapid decline in quasar num- 
bers at low redshift. 
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